UNCLASSIFIED 


Defense  Technical  Information  Center 
Compilation  Part  Notice 

ADP023061 

TITLE:  Radiation  by  a Linear  Array  of  Half-Width  Leaky-Wave  Antennas 
DISTRIBUTION:  Approved  for  public  release,  distribution  unlimited 


This  paper  is  part  of  the  following  report: 

TITLE:  Applied  Computational  Electromagnetics  Society  Journal.  Volume 
21.  Number  3.  November  2006.  Special  Issue  on  Phased  and  Adaptive 
Array  Antennas 

To  order  the  complete  compilation  report,  use:  ADA459097 

The  component  part  is  provided  here  to  allow  users  access  to  individually  authored  sections 
of  proceedings,  annals,  symposia,  etc.  However,  the  component  should  be  considered  within 
the  context  of  the  overall  compilation  report  and  not  as  a stand-alone  technical  report. 

The  following  component  part  numbers  comprise  the  compilation  report: 

ADP023055  thru  ADP023074 


UNCLASSIFIED 


248 


ACES  JOURNAL,  VOL.  21,  NO.  3,  NOVEMBER  2006 


Radiation  by  a Linear  Array  of  Half-Width  Leaky- Wave  Antennas 

(. Invited  Paper) 

Daniel  Killips1,  Joshua  Radcliffe2,  Leo  Kempel1,  and  Stephen  Schneider2 

1 Department  of  Electrical  and  Computer  Engineering,  Michigan  State  University,  East  Lansing, 

MI  48824-1226  (Corresponding  author:  kempel@msu.edu). 

2 Sensors  Directorate,  Air  Force  Research  Laboratory,  Wright-Patterson  AFB,  OH  45433 


Abstract — Leaky-wave  antennas  are  interesting 

apertures  for  a variety  of  applications  due  to  their  low 
profile  and  wide  bandwidth.  They  are  inherently  traveling 
wave  antennas,  and  hence  are  best  suited  for  end-fire 
applications.  A new  type  of  leaky-wave  antenna,  the  half- 
width leaky-wave  antenna  (HWLW),  has  been  recently 
investigated  and  found  to  have  similar  radiation 
properties  as  its  full-width  leaky-wave  (FWLW) 
counterpart,  but  only  requiring  half  the  transverse 
dimension.  In  addition,  the  feeding  mechanism  for  a 
HWLW  antenna  is  considerably  simplified  compared  to 
the  FWLW  antenna.  This  paper  discusses  arraying  these 
antennas  to  provide  both  increased  gain  and  scanning 
capability.  It  will  be  seen  that  arraying  HWLW  antennas 
is  more  complex  than  its  narrowband  counterpart,  the 
patch  antenna. 

Index  Terms — microstrip,  leaky- wave,  transverse 
resonance,  traveling  wave  antenna,  wide  bandwidth, 
array. 

L INTRODUCTION 

MICROSTRIP  leaky-wave  antennas  have  both  low 
profile  and  wide  bandwidth  features.  The  thickness 
of  a microstrip  leaky-wave  antenna  is  no  greater 
than  that  of  the  common  patch  antenna  albeit  with 
substantially  more  bandwidth.  Leaky-wave  antennas  have 
been  extensively  studied  by  Oliner  [1],  Menzel  [2],  Lee 
[3],  and  Lin  [4-6]  among  others.  The  radiating  mode  for 
all  of  these  cases  is  the  EHj  mode,  e.g.  the  first  higher- 
order  mode.  One  of  the  major  issues  with  operating  a 
microstrip  antenna  in  the  first  leaky-wave  mode  involves 
preferential  excitation  of  that  mode  (as  compared  to  the 
fundamental,  non-radiating  EH0  mode).  Menzel  [2] 

utilized  a periodic  array  of  slots  to  suppress  the 
fundamental  mode  while  Lin  et  al.  [4]  utilized  a more 
complex  two-port  feeding  structure.  All  of  the  antennas 
examined  in  these  previous  papers  were  full-width  leaky- 
wave  (FWLW)  microstrip  antennas. 

An  alternative  is  the  half-width  leaky-wave  (HWLW) 
antenna  extensively  studied  by  Thiele  et  al.  [7-9].  This 
antenna  is  formed  by  placing  a shorting  ridge,  between 
the  microstrip  and  ground,  along  one  long  edge  of  the 


antenna.  The  other  long  edge  will  support  a magnetic 
current  wall  that  is  responsible  for  radiation.  Full-width 
microstrip  leaky-wave  antenna  arrays  have  been 
extensively  studied  [10-11].  For  HWLW  antennas,  due  to 
the  relatively  strong  coupling  between  the  array  elements 
via  the  leaky-wave  mode,  it  is  useful  to  utilize  full-wave 
computational  electromagnetics  (CEM)  solvers  rather 
than  to  rely  on  traditional  antenna  array  theory  that 
neglects  coupling.  Since  the  HWLW  antenna  occupies 
half  the  surface  area  of  a conventional  FWLW  antenna, 
the  potential  for  closer  spacing  between  antenna  elements 
is  attractive  for  array  design. 

The  HWLW  antenna  by  itself  has  a major  advantage  over 
the  FWLW  antenna,  relatively  simple  feed  and  loading. 
One  of  the  major  challenges  of  the  FWLW  antenna  is 
suppression  of  the  EH0  mode.  Clearly,  if  the  antenna  is 
operated  at  a frequency  above  the  cut-off  frequency  for 
both  the  EHo  and  EH!  modes,  then  both  modes  can  be 
excited.  Previous  papers  have  discussed  various  methods 
for  suppressing  the  lower-order  mode  in  preference  to  the 
leaky-wave  mode  [2,  4-6,  11].  For  the  HWLW  antenna, 
due  to  the  presence  of  the  shorting  wall  along  one  long 
edge,  the  lower  order  mode  is  automatically  suppressed. 
The  driving  point  impedance  can  be  approximated  using 
an  open  waveguide  model  [6]  along  with  estimated 
wavenumbers  using  the  Transverse  Resonance  Method 
[3].  Using  this  information,  both  the  feed  location  and  the 
location  of  a lumped  load  can  be  estimated.  The  purpose 
of  the  lumped  load  is  to  suppress  the  backward  traveling 
wave  that  will  exist  if  no  load  is  used. 

In  this  present  paper,  the  radiation  properties  of  an  array 
of  HWLW  antennas  are  examined  using  a finite  element- 
boundary integral  (FE-BI)  model.  The  antenna  feeds, 
loads,  and  geometrical  features  will  be  modeled  within 
the  context  of  the  FE-BI  model.  In  addition,  since  this  is  a 
full-wave  model,  the  coupling  between  elements 
attributed  to  the  leaky-wave  mode  will  be  included. 
Radiation  properties  will  be  examined  using  the  model. 

MICROSTRIP  LEAKY-WAVE  ANTENNAS 

A . Leaky-wave  A ntenna  Theory 

The  fundamental  mode  of  a microstrip  line,  the  so-called 
EHo  mode,  is  not  a radiating  mode  (hence  the  popularity 
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of  microstrip  transmission  lines).  The  electric  and 
magnetic  fields  for  this  mode  are  shown  in  Figure  1. 
Radiation  from  such  a structure  can  be  represented  by  two 
long  magnetic  current  walls  (for  the  FWLW  antenna) 
separated  the  width  of  the  microstrip.  Since  these  current 
sources  are  out-of-phase,  the  co-polarized  radiated  field 
in  a plane  bisecting  the  width  of  the  antenna  is  identically 
zero.  Therefore,  as  an  end- fire  antenna,  such  a microstrip 
is  unsuitable. 


Figure  1.  Field  diagram  for  the  EH o mode  (E-field  = 
solid,  H-field  = dashed). 


Rather,  higher-order  modes  must  be  preferentially  excited 
to  realize  a radiating  traveling  wave  structure.  The  first 
higher-order  mode,  the  EHj  mode,  is  one  such  radiating 
mode.  This  mode  (shown  in  Figure  2)  exhibits  electric 
field  odd  symmetry  about  the  axial  centerline  of  the 


Figure  2.  E-field  diagram  for  the  EH!  mode. 


antenna  as  compared  to  the  even  symmetry  of  the  EH0 

mode.  The  radiating  magnetic  currents  are  now  in-phase 
and  hence  radiate  along  the  axis  of  the  antenna. 


The  current  on  either  wall  can  be  represented  as 

i^£)Vz  _/£fe£) 

M(x,y,z)  = ±zA±e  ^ k°  ' e ^ k° 


knz 


(i) 


where  the  attenuation  term  a(z,f)  and  the  propagation 
term  p(z,f)  are  in  general  a function  of  both  position  and 


frequency,  k0  is  the  free-space  wavenumber,  and  the 
wave  coefficients  (a±)  are  associated  with  the  two 
magnetic  wall  currents  at  x = ±w/2  where  wis  the 

width  of  the  full-width  microstrip.  The  attenuation  and 
propagation  terms  for  an  axially  invariant  structure  can  be 
determined  using  the  Transverse  Resonance  Method 
(TRM)  [3].  Once  the  propagation  parameters  are  known, 
the  driving  point  impedance  for  a semi-infinite  line  can  be 


estimated  using  an  open  waveguide  model,  viz.  [5-6]1  as 
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The  input  impedance  for  a finite,  loaded  leaky-wave 
antenna  can  be  estimated  using  the  impedance 
transformation 


= Z„ 


ZL  + zw  tanh(yL) 
ZW+ZL  tanh(yL) 


(3) 


where  ZL  is  the  load  impedance  (for  this  paper, 
ZL  =50Q)and  y = jk  = j(p-ja). 


In  (2),  the  effective  microstrip  width  is  given  by 
Wheeler’s  approximation  [12] 
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Use  of  (2)-(4),  with  the  propagation  parameters  provided 
by  TRM,  allows  determination  of  the  appropriate  feed 
and  load  locations  along  the  width  of  the  strip. 


Assuming  a single  radiating  current  of  the  type 
represented  in  (1),  hence  a leaky-wave  antenna  that  is 
perfectly  terminated  but  aligned  parallel  to  the  x-axis  and 
centered  on  the  origin,  the  radiated  electric  field  in  the 
far-zone  (with  the  spherical  wave  assumed  and 
suppressed)  is  given  by 


E(e,<t>) 


4nL 


[§sin(|) 


+ (|)COS0COS 


(5) 

where  x = jko[sin0cos(j)-k/ko]  and  I0  is  the 


strength  of  the  current.  If  a linear  array  of  uniformly 
spaced,  with  a separation  of  d,  similar  elements  is  aligned 
along  the  y-axis,  the  array  factor  will  be  given  by 

M 


AF=  £. 


>jko(md)sm0sm(|) 


(6) 


m=M 

Obviously,  use  of  (5)  and  (6)  assumes  no  coupling 
between  array  elements. 


B.  Example : Leaky-wave  Antenna  on  Duroid 
As  an  example,  consider  a leaky-wave  antenna  printed  on 
Duroid  5870  (31  mils  thick, er  = 2.33  ,tan5  = 0.0005).  The 
full-width  strip  width  is  15  mm  while  the  half-width  strip 
width  is  7.5  mm.  The  strip  is  taken  to  be  190  mm  long  for 
the  simulations  and  measurements  and  it  is  infinite  for  the 


3 There  is  a typographical  error  in  the  expression  in  [6]  where  the  sine 
function  should  be  squared  as  shown  in  (2).  The  expression  is  correct  in 

[5]. 
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TRM  analysis.  The  propagation  terms,  as  determined 
using  the  transverse  resonance  method,  are  shown  in 
Figure  3.  Traditionally,  the  leaky-wave  region  of 
operation  is  defined  between  the  frequency  such  that 
a - P to  the  frequency  such  that  p = k0 ; hence,  for  the 

example  presented,  approximately  from  6 GHz  to  8 GHz. 
Using  this  information,  along  with  (2)  and  (3),  the  driving 
point  impedance  for  a 190  mm  realization  of  the  antenna 
is  shown  in  Figure  4 where  the  feed  point  is  taken  at  the 
transverse  midpoint  of  the  HWLW  antenna  and  a 50Q 
lumped  load  is  also  placed  along  the  midline  of  the 
antenna  at  the  opposite  end  from  the  feed. 


Frequency  (GHz)) 


Figure  3.  Normalized  attenuation  and  phase  constant  for  a 
leaky-wave  antenna. 


Input  Impedance  for  the  Leaky-Wave  Antenna  [center-fed] 


Figure  4.  Input  impedance  at  the  center  of  the  half-width 
strip  width  for  this  leaky-wave  antenna. 

As  can  be  seen,  the  impedance  seen  at  the  input  port  of 
the  antenna  is  dispersive  and  hence  a perfect  match  to  a 
50  Q line  is  not  possible  across  the  entire  operational 
bandwidth.  In  addition,  as  noted  in  [10,11],  the  traveling 
wave  fields  guided  by  a microstrip  structure  couple  rather 


strongly  to  the  leaky-waves  emanating  from  neighboring 
microstrips  causing  deviation  from  theoretical  predictions 
for  both  the  wavenumber,  driving  point  impedance,  and 
consequently,  the  radiated  field.  Hence,  use  of  a full-wave 
computational  method  is  desirable  for  designing  leaky- 
wave  arrays. 

II.  FINITE  ELEMENT-BOUNDARY  INTEGRAL 
MODEL 

The  hybrid  finite  element-boundary  integral  (FE-BI) 
method  is  widely  understood  to  be  an  accurate  and 
efficient  method  for  modeling  planar,  microstrip 
apertures.  The  assumed  geometry  is  a cavity-backed 
aperture  recessed  within  an  infinite  metallic  plane.  The 
FE-BI  formulation  is  an  integro-differential  equation 
representing  a weak  enforcement  of  the  relevant  boundary 
conditions  applied  to  the  vector  wave  equation.  This  FE- 
BI  expression  is  given  as 

[ [(V  X W( ) • pr-’ ■ (v  xJE-nt )- W,  ■ fr  • ] dV  + 

k02  ||Wj  • [z  x <fe2  x zj-  E”11  dS'dS  = 

- jk0Z0  £ Wj  • J““p  dV 

(7) 

where  W,  is  a vector  testing  function,  Ge0  is  the  half- 
space dyadic  Green’s  function,  and  Jimpis  an  impressed 
current  used  to  represent  the  antenna  feed.  The  material  in 
the  computational  volume,  e.g.  the  antenna  cavity,  is 
represented  by  the  relative  permittivity  tensor  (er)  and 
the  relative  inverse  permeability  tensor (jl"1 ).  The  FE-BI 

formulation  is  reduced  to  a linear  system  via  Galerkin’s 
procedure  resulting  in  the  following  linear  system 


In  this,  the  unknown  coefficients  associated  with  the 
expansion  functions  (Wj ) are  represented  by  )-  For  the 

unknown  coefficients,  (e^)  represents  unknowns 

associated  with  the  surface  while  {eJ  ) represent 

unknowns  not  associated  with  the  surface.  This  system 
can  efficiently  be  solved,  both  in  terms  of  memory 
consumption  and  wallclock  time,  using  an  iterative  solver 
such  as  the  biconjugate  gradient  method  [14]. 

By  utilizing  a fast  Fourier  transform  (FFT)  method  to 
accelerate  the  computation  of  matrix-vector  products  in 
an  iterative  solution  of  the  FE-BI  linear  system, 
simulation  of  a large  antenna  is  feasible  within  practical 
design  time  limits.  Hence,  brick  elements  were  used  to 
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discretize  the  computational  volume  [13].  The  shorted 
side  of  the  HWLW  antenna  is  modeled  using  an 
infinitesimally  thin  perfectly  conducting  wall  running 
from  one  end  of  the  strip  to  the  other  end.  The  lumped 
loads  were  modeled  as  infinitesimally  thin  loads  placed 
along  the  shorted  edge  of  the  half-width  antenna.  The 
feeds  are  modeled  likewise  as  infinitesimally  thin  probe 
feeds  placed  along  the  shorted  edge  of  the  antenna 
opposite  the  loads  on  the  mid-line  of  the  antenna. 

In  addition,  the  computer  program  was  parallelized  using 
OpenMP  [15]  and  run  on  a dual-core  AMD  64-bit  system 
with  a total  of  four  cores  per  node.  The  program  has  also 
been  used  on  an  SGI  Altix  3700  Bx2  using  64-bit  Itanium 
processors  [16].  The  primary  routine  that  was  parallelized 
was  the  iterative  matrix  solvers  (the  biconjugate  gradient 
and  transpose-free  QMR  algorithms  are  implemented) 
and  fast  Fourier  transforms.  For  this  work,  the  BiCG 
algorithm  [14]  was  used. 


HI.  COMPUTED  DATA  FOR  A FIVE  ELEMENT 
LINEAR  ARRAY 

To  validate  the  computational  model  discussed  above,  an 
experiment  was  conducted  in  a compact  anechoic 
chamber  using  a single  HWLW  antenna.  The  shorting 
wall  in  the  half-width  antenna  was  approximated  using 
shorting  pins  that  were  5 mm  apart  running  the  length  of 
the  antenna  strip.  The  ports  were  attached  to  an  HP-8510 
vector  network  analyzer.  A representative  radiation 
pattern,  comparing  computed  and  measured  data,  is 
shown  in  Figure  6. 


Normalized  Gain  [6.8  GHz] 


Figure  6.  Comparison  of  computed  and  measured 
radiation  patterns  from  a half-width,  leaky-wave  antenna 
at  6.8  GHz. 


The  agreement  is  reasonably  good  considering  that  the 
shorting  wall  was  not  solid,  the  FE-BI  model  used  a probe 
feed,  and  the  actual  groundplane  was  finite  while  the 
modeled  one  was  infinite. 

Next,  the  FE-BI  model  was  used  to  determine  the  impact 
of  leaky-wave  coupling  on  the  fields  supported  by 
adjacent  microstrip  elements.  Figure  7 illustrates  a 
representative  coupling  between  a driven  element  (the 
lower  one  in  the  figure)  and  an  unfed  element  that  is 
terminated  on  both  ends  of  the  antenna  with  50D  loads 
and  1.95  cm  (center-to-center)  away  from  the  lower 
element.  The  simulation  was  conducted  with  a single 
cavity  containing  both  the  driven  element  and  the  passive 
element.  In  this  manner,  coupling  between  elements  can 
occur  via  substrate,  across  the  surface  of  the  aperture,  and 
via  space-wave  mechanisms.  For  an  aperture  where  each 
element  is  contained  within  individual  cavities,  coupling 
via  the  substrate  is  suppressed.  Since  this  is  assumed  to  be 
the  dominant  mechanism,  since  the  EHj  mode  is  a leaky- 
wave  mode,  this  geometrical  arrangement  is  appropriate. 


-10 


-10  -5  0 5 10  HI  4q 

Figure  7.  Normal  fields  in  cavity  for  driven  lower  element 
and  passive  upper  element. 

It  is  clear  that  the  leaky-wave  from  the  driven  lower 
element  is  coupled  into  the  upper  element.  Notice  that  the 
fields  in  the  upper  element  are  strongest  in  the  end  of  the 
element  opposite  the  driving-point  in  the  lower  element;  a 
phenomena  expected  from  leaky-wave  coupling  since  the 
driven  leaky-wave  will  have  a directional  component 
from  left  to  right  in  the  figure.  To  assess  the  impact  of 
this  coupled  field  to  the  radiating  current  (and  hence  the 
fidelity  of  the  radiation  pattern  to  predictions  based  on  a 
single  traveling  wave),  the  coupling  ratio  was  determined 
for  various  separations.  The  coupling  ratio  is  defined  as 
the  ratio  of  the  maximum  field  strength  in  the  region  of 
the  passive  antenna  to  the  maximum  field  strength  in  the 
region  of  the  driven  element.  Hence,  it  is  a measure  of 
how  strong  one  leaky-wave  couples,  and  hence  corrupts, 
the  leaky- wave  propagating  in  another  microstrip.  As  a 
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function  of  the  center-to-center  separation  distance,  the 
coupling  between  these  elements  at  7 GHz  is  given  in  the 
Figure  8.  As  can  be  seen,  the  coupling  remains  fairly 
strong  even  out  to  a distance  of  6.45  cm.  For  7 GHz,  the 
free-space  wavelength  is  4.29  cm  and  hence  a separation 
of  6.45  cm  would  result  in  grating  lobes.  This  result 
enforces  the  assumption  made  previously  that  to 
accurately  assess  the  radiation  properties  of  a leaky-wave 
array,  either  the  rather  complex  formulation  presented  in 
[1 1]  or  a full-wave  model  such  as  the  FE-BI  method  must 
be  used. 


Figure  8.  Coupling  ratio  between  elements  vs.  separation 
distance. 


A five-element  array  was  simulated  using  the  FE-BI 
model  and  having  a separation  of  1 .2  cm  between 
elements  (e.g.  less  than  1/3 rd  of  the  free-space  wavelength 
to  avoid  grating  lobes).  The  theoretically  predicted 
radiation  patterns  based  on  (5)  and  (6)  - e.g.  the 
normalized  radiation  intensity  [17]  - is  shown  in  Figures 
9 and  10  for  a uniformly  illuminated  array. 


Radiation  Pattern  for  ^component 


4>[deg] 

Figure  9.  Theoretical  (J)-(J)  radiation  pattern  for  a uniformly 
illuminated  array. 


..Radiation  Pattern  for  80-component 
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Figure  10.  Theoretical  0-0  radiation  pattern  for  a 
uniformly  illuminated  array. 

The  corresponding  patterns  as  computed  by  the  FE-BI 
model  are  shown  in  Figures  1 1 and  12,  respectively. 


Radiation  Pattern  for  ^component 
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Figure  11.  Numerical  radiation  pattern  for  a 
uniformly  illuminated  array. 

For  all  of  the  contour  plots,  each  contour  represents  3 dB 
of  change.  As  seen  by  comparing  Figure  9 to  Figure  11 
and  Figure  10  to  Figure  12,  the  agreement  is  quite  good 
with  the  exception  of  the  somewhat  diffused  sidelobes  for 
the  numerically  computed  results.  The  reason  is 
presumably  due  to  the  leaky-wave  coupling  that  is 
neglected  in  the  theoretical  results. 

Another  example  is  the  same  array;  however,  the  main 
beam  is  now  steered  30  degrees  from  the  array  axis 
normal.  The  theoretical  results  are  shown  in  Figures  13 
and  14  and  the  numerical  results  are  shown  in  Figures  15 
and  16,  respectively. 
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Figure  12.  Numerical  0-0  radiation  pattern  for  a 
uniformly  illuminated  array. 


Figure  13.  Theoretical  <J>-(J)  radiation  pattern  for  an  array 
steered  30  degrees  off  the  x-axis. 


Figure  14.  Theoretical  0-0  radiation  pattern  for  an  array 
steered  30  degrees  off  the  x-axis. 


Figure  15.  Numerical  radiation  pattern  for  an  array 
steered  30  degrees  off  the  x-axis. 
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Figure  16.  Numerical  0-0  radiation  pattern  for  an  array 
steered  30  degrees  off  the  x-axis. 

Again,  the  agreement  between  the  theoretical  and 
numerical  results  is  quite  good  with  the  exception  of 
higher  sidelobes  for  the  numerical  results. 


IV.  Conclusions 

Microstrip  leaky-wave  antennas  offer  an  attractive  end- 
fire  aperture  that  is  thin  and  easy  to  manufacture.  The 
half-width  leaky-wave  antennas  investigated  herein  have 
a particularly  simplified  feed  structure  compared  to  the 
more  traditional  full-width  leaky- wave  antennas.  An  array 
of  these  antennas  was  investigated  using  a hybrid  finite 
element-boundary  integral  method.  Strong  coupling 
between  adjacent  elements  was  observed  and  quantified. 
Nevertheless,  the  radiation  patterns  between  the 
theoretical  and  numerical  results  agreed  quite  well  except 
for  the  sidelobes.  It  is  assumed  that  the  leaky-wave 
coupling  between  elements  perturbs  the  radiating 
magnetic  current  for  each  element  causing  changes  in  the 
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sidelobes.  Hence,  for  applications  requiring  tight  sidelobe 
control,  use  of  a full-wave  analysis  method  is  essential. 
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